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Abstract 

Using full-potential density functional calculations we have investigated the structural 
and electronic properties of graphene and some of its structural analogues, viz., mono- 
layer (ML) of SiC, GeC, BN, A1N, GaN, ZnO, ZnS and ZnSe. While our calculations 
corroborate some of the reported results based on different methods, our results on ZnSe, 
the two dimensional bulk modulus of ML-GeC, ML-A1N, ML-GaN, ML-ZnO and ML-ZnS 
and the effective masses of the charge carriers in these binary mono-layers are something 
new. With the current progress in synthesis techniques, some of these new materials may 
be synthesized in near future for applications in nano-devices. 
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1 Introduction 



Graphene is a crystal of carbon (C) atoms tightly bound in a two dimensional (2D) hexag- 
onal lattice. It is a monolayer of carbon atoms (ML-C), i.e., one atom thick. The exotic 
properties of this 2D material were revealed only in 2004-2005 by a series of papers com- 
ing from the Manchester [IH3] and Columbia [I] groups. The unambiguous synthesis (by 
mechanical exfoliation of graphite), identification (by transmission electron microscopy 
(TEM)) and experimental determination of some of the exotic properties of graphene were 
reported first in 2004, by the Manchester group led by Novoselov and Geim pQ. In 2010, 
Konstantin S. Novoselov and Andre Geim were awarded the Nobel Prize in Physics for the 
"groundbreaking experiments regarding the two-dimensional material graphene". How- 
ever, graphene research has a history which dates back to the 1859 work of Brodie |5HH] - 
The term graphene was introduced by Boehm and his colleagues in 1986 [9] . Graphene = 
"graph" + "ene" and the term "graph" is derived from the word "graphite" and the suf- 
fix "ene" refers to polycyclic aromatic hydrocarbons. Now graphene is considered as the 
basic building block of graphitic materials (i.e., graphite = stacked graphene, fullerenes 
= wrapped graphene, nanotube = rolled graphene, graphene nanoribbon = nano-scale 
finite area sized rectangular graphene). 

The existence of free-standing 2D crystals were believed impossible for several years, 
because they would ultimately turn into a three-dimensional (3D) objects as predicted 
by Peierls [10], Landau [TT] and Mermin [12]. The theoretical reason for this is that 
at finite temperatures, the displacements due to thermal fluctuations could be of the 
same magnitude as the inter-atomic distances, which make the crystal unstable. Fur- 
ther, experimentally one generally finds that thin films cannot be synthesized below a 
certain thickness, due to islands formations or even decomposition. Hence, the synthesis 
of graphene [1] was surprising which put a question mark on the predictions of Peierls, 
Landau and Mermin. However, this issue was (at least partially) solved when it was shown 
that freestanding graphene sheets display spontaneous ripples owing to thermal fluctua- 
tion [US], and therefore real graphene is not perfectly flat. It is important to note that 
such instabilities are the result of thermal fluctuations which disappear at temperature T 
= 0° K. This aspect will be used later in our discussions on the study of the stability of 
graphene and some of its structural analogues reported here on the basis of their ground 
state (i.e., T = 0° K) total energies. 

Being a one-atom-thick planar crystal of C atoms, graphene is the thinnest nano- 
material ever known. It has exotic mechanical, thermal, electronic, optical and chemical 
properties, such as the high carrier mobility, a weak dependence of mobility on carrier con- 
centration and temperature, unusual quantum hall effect, hardness exceeding 100 times 
that of the strongest steel of same thickness and yet flexible (graphene can sustain elastic 
tensile strain more than 20% without breaking, and is brittle at certain strain limit), high 
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thermal conductivity comparable to that of diamond and 10 times greater than that of 
copper, negative coefficient of thermal expansion over a wide range of temperature. Be- 
cause of these properties graphene has potentials for many novel applications [6l[7l fT4H20] . 
The rapid advancements of nanotechnology and the computing power have enabled the 
researchers in the field to explore the unusual properties of graphene from many per- 
spectives of application and fundamental science. The study of graphene is possibly the 
largest and fastest growing field of research in material science. 

The impressive growth in the research on graphene has inspired the study of other 
graphene-like 2D materials [20|l2~T]. For instance, a number of 2D/quasi-2D nanocrystals 
(not based on carbon) have been synthesized or predicted theoretically in recent years. 
Representative samples of other 2D nanocrystals which have been synthesized include 
BN, MoS 2 , MoSe 2 , Bi 2 Te 3 USEE], Si [22], ZnO [23]. Recently, the density functional 
theory (DFT) calculations of Freeman group [23] have shown that when the layer num- 
ber of (OOOl)-oriented wurtzite (WZ) materials (e.g., A1N, BeO, GaN, SiC, ZnO and 
ZnS) is small, the WZ structures transform into a new form of stable hexagonal BN-like 
structure. This prediction has recently been confirmed in respect of ZnO [23], whose sta- 
bility is attributed to the strong inplane sp 2 hybridized bonds between Zn and O atoms. 
Graphene-like 2D / quasi-2D honeycomb structures of group- and /// — V binary com- 
pounds have also been studied [251126] by using pseudo-potential DFT calculations. Here, 
we report our calculations on the structural and electronic properties of graphene and 
some other graphene-like monolayer (ML) structures of the binary compounds, viz., SiC, 
GeC, BN, A1N, GaN, ZnO, ZnS, ZnSe using the DFT. 

2 Computational Methods 

We use the DFT based full-potential (linearized) augmented plane wave plus local orbital 
(FP-(L)APW+lo) method [271429] as implemented in the elk-code (http: / / elk.sourceforge.net/) 
and the Perdew-Zunger [30] variant of local density approximation (LDA) for our calcu- 
lations. The accuracy of this method and code has been successfully tested in our pre- 
vious studies [3TH3T] . For plane wave expansion in the interstitial region, we have used 
8 < |G + k| max . x R mt < 9, where R mt is the smallest muffin-tin radius, for deciding the 
plane wave cut-off. The Monkhorst-Pack [38J fc-point grid size of 20 x 20 x 1 was used 
for structural and of 30 x 30 x 1 for band structure and total density of states (TDOS) 
calculations. The total energy was converged within 2 /zeV/atom. The 2D-hexagonal 
structure was simulated by a 3D-hexagonal super cell with a large value of c-parameter 
(= |c| =40 a.u.) as shown in Figure 1. For structure optimization, we have considered 
two different structures, viz., (i) Planar structure (PL) and (ii) Buckled structure (BL) 
as shown in Figure 1 for a monolayer of BN (ML-BN), which is the prototype of all the 
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(e) PL - Unit super cell (f) BL - Unit super cell 



Figure 1: Two possible structures of ML-BN: Planar (PL) and Buckled (BL) in ball- 
and-stick model. Sub-figures (a) and (b) show the top-down views of ML-BN in PL 
and BL structure respectively. In the buckled ML-BN, B atoms (large balls, red) and N 
atoms (small balls, green) are in two different parallel planes (sub-figure (d)), buckling 
parameter A is the perpendicular distance between those parallel planes and A = A for 
a flat ML-BN (sub- figure (c)). Sub-figures (e) and (f) respectively depict the 3D super 
cell of hexagonal BN in PL and BL structure with a large value of c-parameter (the length 
along the vertical axis). 

materials considered here (for graphene, B and N atoms are to be replaced by the C 
atoms). 

3 Results and Discussions 

We have optimized the 2D hexagonal structures of graphene (ML-C) and the mono- 
layers of SiC (ML-SiC), GeC (ML-GeC), BN (ML-BN), A1N (ML-A1N), GaN (ML-GaN), 
ZnS (ML-ZnS) using the principle of minimum energy for the stable structure as per 
the following procedure. Initially we assumed the planar structures for all the materials 
considered here. With this assumption, we calculated the ground state in-plane lattice 
parameter (a ) of each of these structures as listed in Table-1 along with the available 
reported values by other authors. Then we investigated the possibility of buckling in 
these structures at their assumed planar ground states (characterized by their respective 
a values in planar states) by introducing the concept of bucking (Figure 1) and the 
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principle of minimum energy for the most stable structure. Our calculated variation of 
total energy (E) with buckling parameter (A) at fixed value of was then plotted to 
look for the value of A which corresponds to the minimum energy. 




t-l ■ 1 1 1 1 1 — 

0.00 0.05 0.10 0.15 

Buckling parameter: A (A) 



Figure 2: Buckling probe of monolayer graphene (ML-C). 

In Figures 2-4, we present our calculated results selecting one example from each 
group of materials considered here: graphene (ML-C) for Group-IV material, ML-BN for 
Group-III-V material and MLZnS for Group-II-VI material. As seen in the Figures 2-4, 
these structures have minimum energy at A = 0.00 A , which means that these materials 
adopt 2D planar structures in the ground state (T = 0° K) unlike the case with silicene 
(graphene analogue of Si) (22l[3Tj that adopts a buckled structure in its ground state. It 
is to be noted that these results do not conflict with the theory of stability of 2D crys- 
tals [10HT2]. Our calculated optimized structural parameters of the materials are listed in 
Table 1 along with available reported results based on other methods. The stability tests 
ML-ZnO and ML-ZnSe are not yet complete. 

Considering the facts that LDA usually underestimates and GGA [39] usually over- 
estimates the lattice constant, the calculated value of ao slightly depends on the c- 
parameter used for super cells in DFT based calculations as we have seen before [31] 
for graphene (0.12% lower value of a for c — > 00) and silicene (0.075% lower value of a 
for c 00), and the different methods of study used by different authors, our results in 
Table 1 are acceptable. 

Using the definition of the 2D bulk modulus of planar a graphene-like 2D material as 
B2D = A(d 2 E I dA 2 )\A min (where A is the area of the periodic cell of the 2D lattice and 
A min is the area with minimum energy), we have calculated the bulk modulus of each 
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Table 1: Calculated LDA values of ground state in-plane lattice constant ao(— l a l — 
|b|), buckling parameter A of monolayer graphene (MLC) and some structurally similar 
binary compounds compared with available calculated data. PAW stands for the projector 
augmented wave. 



Material 


a (A) 


A (A) 


Remark / Reference 




z.44o 


U.U 


Uur Work/ |oi4ljj] 


1\/TT H 


AP. 


n n 

U.U 


r/vvv-poLenxiai/ [zo] 




2.4431 


0.0 


Pseudo-potential/ [26J 




3.066 


0.0 


Our Work/ [37] 


ML-SiC 


3.07 


0.0 


PAW-potential/ [25] 




3.0531 


0.0 


Pseudo-potential/ [26J 


ML-GeC 


3.195 


0.0 


Our Work/ [37] 




3.22 


0.0 


PAW-potential/ [25] 




2.488 


0.0 


Our Work/ [33l|3l] 


ML-BN 


2.51 


0.0 


PAW-potential/ [25] 




2.4870 


0.0 


Pseudo-potential/ [26J 


ML-A1N 


3.09 


0.0 


Our Work/ [31] 




3.09 


0.0 


PAW-potential/ [25] 


ML-GaN 


3.156 


0.0 


Our Work/ [37] 




3.20 


0.0 


PAW-potential/ [25] 


ML-ZnO 


3.20 


0.0 (used) 


Our Work/ [35] 




3.283 


0.0 


PAW-potential/ [46] 


ML-ZnS 


3.7995 


0.0 


Our Work/ [36] 




3.890 


0.0 


PAW-potential with GGA [39]/ [47] 


ML-ZnSe 


3.996 


0.0 (used) 


Our Work/ [37J 
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Table 2: Calculated 2D bulk modulus (B 2 d) value of monolayer graphene and some 
other mono-layers of binary compounds in graphene-like structure compared with available 
calculated data using other methods. Our cited values of B 2 d from Wang [26] were 
obtained from Wang's calculated values of 2D elastic constants 711 and 712 and another 
definition of B 2 d = (711 + 712 )/2 [4"9] . 



Material (2D) 


B 2D (N/m) 


Remark / Reference 


ML-C 


223.85 


This Work 




214.41 


Pseudo-potential/ [26J 




220.00 


Pseudo-potential/ [48J 


ML-SiC 


125.66 


This Work 




121.945 


Pseudo-potential/ [26J 


ML-GeC 


113.39 


This Work 


ML-BN 


188.03 


This Work 




181.91 


Pseudo-potential/ [26J 




212.0 


Born's Perturbation Method/ [49J 


ML-A1N 


113.44 


This Work 


ML-GaN 


109.45 


This Work 


ML-ZnO 


94.78 


This Work 


ML-ZnS 


53.94 


This Work 


ML-ZnSe 


48.57 


This Work 
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0.00 0.05 0.10 0.15 
Buckling parameter: A (A) 



Figure 4: Buckling probe of monolayer ZnS (ML-ZnS). 

of the materials considered here. These calculated values of B2D are listed in Table 2 
along with the available calculated data based on other methods. As seen in Table 2, 
the 2D bulk modulus of graphene has the highest value and that of ML-ZnSe has the 
lowest value among the materials considered here, implying the fact that graphene has 
stronger in-plane bonds than the others and ML-ZnSe has the weakest bond among these 
materials. ML-GeC is seen to have same strength as that of ML-A1N. These data may be 
useful in selecting the 2D materials having desired level of strength in nano-mechanical 
applications in specific situations. 

The LDA band structures of graphene (ML-C) and some other mono-layers of SiC, 
GeC, BN, A1N, GaN, ZnO, ZnS and ZnSe in graphene-like planar structure are depicted 
in Figure 5 and in Table 3, we have listed our calculated LDA values of the band gaps of 
these materials along with the reported values. As seen in the explicit case of ML-BN in 
Table 3, although LDA under-estimates the band gap (in our case the calculated value is 
about 23% less than the experimental value of 5.971 eV) it correctly predicts the nature 
of the gap. However, the use of a more advanced approximation, such as the LDA + GWo 
used in [25], improves the band gap problem (in the case of ML-BN although it overesti- 
mates band gap by 15% in this case), but the authors of [25j have presented their result 
in terms of an indirect band gap for ML-BN as opposed to the experimental observation 
of direct band gap in ref . [10] . In all other cases, the LDA nature of band gap is the same 
as that of the LDA + GWo band gap. The important point is that the actual band gap 
is always more than the LDA value. 
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Table 3: Calculated LDA band gaps compared and listed with reported values. 



Material (2D) 


E 9 (eV) 


Remark / Reference 


ML-C 


(if if) 
(if if ) 



Our Work/ [31,33] 
PAW-potential/ [25] 
Pseudo-potential/ [26J 


ML-SiC 


2.547(if -» M) 
2.52 (if M) 
4.19 (if — )■ M) 
2.39 


Our Work/ [37] 
PAW-potential/ [25] 
PAW-potential (LDA + GW )/ [25] 
Pseudo-potential/ [26J 


IV TT 

ML-GeC 


2.108(if — > A) 
2.09 (if if) 
3.83 (if — )■ a ) 


Our Work/ [37] 
PAW-potential/ [25] 
PAW-potential (LDA + GW )/ [25] 


A /TT "D AT 


4.606(A — > K ) 
4.61 (if -> if) 
6.86 (if -> r) 
4.35 

5.971 (direct) 


Uur Work/ [33,34J 
PAW-potential/ [25] 
PAW-potential (LDA + GW )/ p5] 
Pseudo-potential/ [26J 
Experiment/ [40] 


ML-A1N 


3.037(if -> T) 
3.08 (if -> r) 
5.57 (If -> r) 


Our Work/ [34] 
PAW-potential/ [25] 
PAW-potential (LDA + GW )/ [25] 


ML-GaN 


2.462(if -> r) 
2.27(if -> r) 


Our Work/ [37] 
PAW-potential/ [25] 


ML-ZnO 


i.680(r -> r) 
i.68 (r ->■ r) 


Our Work/ [35] 
PAW-potential/ [46] 


ML-ZnS 


2.622(r -> r) 
2.07 (r -> r) 


Our Work/ [36] 

PAW-potential with GGA [39]/ [47] 


ML-ZnSe 


1.866 (r -> r) 


Our Work/ [37] 
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Table 4: Calculated values of effective masses of the charge carriers in ML-C, ML-Si, 
ML-Ge, ML-BN, ML-A1N, ML-GaN, ZnO, ZnS, ZnSe. m is the rest mass of electron. 



Material (2D) 


Effective mass (m ) 


Remark / Reference 


ML-C 


m e (K) = 


This Work 




m h {K) = 


This Work 


ML-SiC 


m e (M) = 0.411 


This Work 




m h {K) = -0.488 


This Work 


ML-GeC 


m e (K) = 0.509 


This Work 




m h {K) = -0.400 


This Work 


ML-BN 


m e (K) = 0.920 


This Work 




m h (K) = -0.617 


This Work 


ML-A1N 


m e (T) = 0.523 


This Work 




m h {K) = -1.470 


This Work 


ML-GaN 


m e (r) = 0.266 


This Work 




m h (K) = -1.157 


This Work 


ML-ZnO 


m e (r) = 0.253 


This Work 




m lh {Y) = -0.374 


This Work 




m hh {Y) = -0.793 


This Work 


ML-ZnS 


m e (T) = 0.173 


This Work 




m lh {Y) = -0.154 


This Work 




mw,(r) = -0.665 


This Work 


ML-ZnSe 


m e (T) = 0.107 


This Work 




m, h (r) = -0.103 


This Work 




m hh (r) = -0.652 


This Work 
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Figure 5: Band structures of graphene (ML-C) and some other mono-layers of SiC, GeC, 
BN, A1N, GaN, ZnO, ZnS and ZnSe in graphene-like planar structure within LDA. 



In Table 4, we present our calculated values of the effective masses of the charge car- 
riers in ML-C, ML-Si, ML-Ge, ML-BN, ML-AIN, ML-GaN, ZnO, ZnS, ZnSe, determined 
at the band edges at the special points as appropriate for the material under study. From 
the linear energy dispersion close to the K point of the hexagonal Brillouin Zone (BZ) as 
shown in Figure 6, mass-less carriers in graphene were inferred. As seen in Figure 6, near 
the K point of the BZ, the energy bad dispersion is linear in k: 

E± = ±v F hk (1) 

where vp is the magnitude of Fermi velocity of the charge carriers in graphene and (hk) is 
the magnitude of the momentum. Thus, charge carriers in graphene behave like mass-less 
relativistic particles. From the slope of the linear bands one obtains the value of vf, which 
in our calculation corresponds to vp = 0.8327 x 10 6 m/s. This value of Vp is close to the 
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Figure 6: Linear energy dispersion near the K point in monolayer graphenes band struc- 
ture. 

experimentally measured value of vp = 0.79 x 10 6 m/s in graphene monolayer deposited on 
graphite substrate |JT] and the reported calculated LAPW value of 0.833 x 10 6 m/s [12] 
and about 20% less than that measured in coupled multi-layers [HJIH] and the tight- 
binding value: Vp = 1 x 10 6 m/s [TS]. The smallness in the value of vp in ref. [H] was 
attributed to the electron-phonon interactions due to strong coupling with the graphite 
substrate. However, the closeness of our calculated value of vf for a freestanding graphene 
monolayer with the results of ref. [H] suggests a very weak coupling of graphene to the 
graphite substrate used in ref. [U]. Surprisingly, our calculated value of Vp is close to the 
measured value of Vp = 0.81 x 10 6 m/s [15] in metallic single-walled carbon nanotube. 

4 Conclusions 

Using full potential DFT calculations we have investigated the structural and electronic 
properties of graphene and some other graphene-like materials. While our results cor- 
roborate the previous theoretical studies based on different methods, our calculations on 
ML-ZnSe, the two-dimensional bulk modulus of ML-GeC, ML-A1N, ML-GaN, MLZnO, 
ML-ZnS, and the effective masses of the binary monolayer compounds considered here 
are our new results. We hope, with the advancement of fabrication techniques, the hypo- 
thetical graphenelike materials discussed here will be synthesized in the near future for 
potential applications in a variety of novel nano-devices. 
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